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Abstract 

We study the local volatility function in the Foreign Exchange market where both domestic and 
foreign interest rates are stochastic. This model is suitable to price long-dated FX derivatives. We 
derive the local volatility function and obtain several results that can be used for the calibration 
of this local volatility on the FX option's market. Then, we study an extension to obtain a more 
general volatility model and propose a calibration method for the local volatility associated to this 
model. 

Keywords : Local volatility; Stochastic volatility; Foreign Exchange; Stochastic interest rates; 
Calibration. 



1 Introduction 



Recent years, the long-dated FX option's market has grown considerably. Currently most traded 
and liquid long-dated FX Hybrid products are Power-Reverse Dual-Currency swaps (PRDC) (see for 



example Piterbarg, 2006| ) as well as vanilla or exotic long-dated products such as barrier options. 
While for short-dated options (less than 1 year), assuming constant interest rates does not lead to 
significant mispricing, for long-dated options the effect of interest rate volatility becomes increasingly 
pronounced with increasing maturity and can become as important as that of the FX spot volatility. 



Most of the dealers are using a three-factor pricing model for long-dated FX products (see Piterbarg 



2006, Sippel and Ohkoshi 2002| ) where the FX spot is locally governed by a geometric Brownian 
motion, while each of the domestic and foreign interest rates follows a Hull- White one factor Gaussian 
model |Hull and White" , 1993 . Using such a model does not allow the volatility smile/skew effect 
encountered in the FX market to be taken into account, and is therefore not appropriate to price and 
hedge long-dated FX products. 



Different methods exist to incorporate smile/skew effects in the three-factor pricing model. In 
the literature, one can find different approaches which consist either of using a local volatility for the 
FX spot or a stochastic volatility and/or jump. There are many processes that can be used for the 
stochastic volatility and their choices will generally depend on their tractability and solvability. All 



these models should be calibrated over the market before being used for pricing. The calibration is 
normally based on calculating prices of liquid products for different strikes and maturity and the pa- 
rameters of the model are adjusted until these prices match sufficiently with the market. However, in 
most cases it is difficult to derive analytical formulae, and consequently the calibration procedure often 
remains approximative or computationally demanding. Andreasen suggested in Andreasen 2006 1 a 



stochastic volatility approach by combining a Heston |Heston 1993 stochastic volatility model with 
independent stochastic interest rates. He has derived closed-form Fourier expressions for vanilla op- 
tions which are useful for calibration on the vanilla market smile/skew. He has used an indirect 
approach in the form of a volatility displacement parameter to correlate the independent interest 
rates with the spot FX rate. In |Antonov et al. 2008 , Antonov et al. have underlined the problem 
that using non-trivial correlations destroys the affine structure and exact solvability. By using the 
technique of Markovian Projection, they have derived approximation formulae for the calibration of 
FX options in a three-factor model coupled with a Heston stochastic volatility under a full correlation 
structure. In Ivan Haastrecht et al. 2009 , van Haastrecht et al. have obtained the exact pricing of 



FX options under the three-factor model coupled with a Schobel and Zhu Schobel and Zhu, 1999 



stochastic volatility and the full correlation structure. Their model could cover Poisson type jump with 



a trivial extension. Recently, van Haastrecht and Pelsscr considered in van Haastrecht and Pelsser 



2009] the pricing of FX options under Schobel and Zhu and also the Heston stochastic volatility with 
a multi- factor Gaussian interest rates and a full correlation structure. More precisely, they derived 
the characteristic functions required for the Fourier-based pricing methods under Schobel and Zhu 
stochastic volatility. Unfortunately, they did not obtain a closed-form expression for the characteris- 
tic functions under the Heston stochastic volatility with a full correlation structure. However, they 
presented a calibration method based on the characteristic functions in the special uncorrelated case 
by using a projection of the general model onto the uncorrelated case, or by using it as a control 
variate for the general model. Using martingale methods and Fourier inversion techniques, Ahlip 
Ahlip, 20081 has derived an analytical formula for the price of a European call on the spot FX rate 



in the case of a stochastic volatility model with stochastic interest rates where domestic and foreign 
interest rates are modelled by Ornstein-Uhlenbeck processes and the instantaneous volatility follows 
a mean-reverting Ornstein-Uhlenbeck process correlated with the spot FX rate. Finally, L.A. Grzelak 



and C.W. Oosterlee Grzelak and Oosterlee, 2010 have derived semi-closed form approximations for 



the forward characteristic function in a foreign exchange model of Heston-type, in which the domestic 
and foreign interest rates are generated by the short-rate process of Hull- White and have extended the 
framework by modeling the interest rate by a stochastic volatility displaced-diffusion Libor Market 
Model. 



A local volatility framework was explored by Piterbarg in |Piterbarg, 2006 where the volatility of 



the spot FX rate was a function of both the time and the spot itself. He has derived an approxima- 
tive formula for the local volatility which allows for a fast calibration of the model on vanillas. The 
calibration essentially captures the "slope" of the implied volatility surface but does not exactly fit its 
convexity. 



Local volatility models, introduced in 1994 by Dupire |Dupire 1994] and Derman and Kani |Der- 
man and Kani , 19941 in the settings of equity, have the benefit over a stochastic volatility model that 



they are Markovian in only one factor since the local volatility is a deterministic function of both 
the FX spot and time. It avoids the problem of working in incomplete markets in comparison with 
stochastic volatility models and is therefore more appropriate for hedging strategies. Local volatil- 
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ity models also have the advantage to be calibrated on the complete implied volatility surface, and 
consequently local volatility models usually capture more precisely the surface of implied volatilities 
than stochastic volatility models. However, a local volatility model has the drawback that it predicts 
unrealistic dynamics for the stock volatility since the volatilities observed in the market are really 
stochastic, capable of rising without a movement in spot FX prices. In |Bossens et al. 2010 , the 
authors compare short-dated barrier option market prices with the corresponding prices derived from 
either a Dupire local volatility or a Heston stochastic volatility model both calibrated on the vanilla 
smile/skew. It appears from that study that in a simplified world where exotic option prices are 
derived either from Dupire local volatility or from Heston stochastic volatility dynamics, a FX market 
characterized by a mild skew (USDCHF) exhibits mainly a stochastic volatility behavior, and that 
FX markets characterized by a dominantly skewed implied volatility (USDJPY) exhibit a stronger 
local volatility component. This observation also underlines that calibrating a stochastic model to 



the vanilla market is by no means a guarantee that exotic options will be priced correctly Schoutens 
et al. , 2004 , as the vanilla market carries no information about the smile dynamics. The market 
dynamics could be better approximated by a hybrid volatility model that contains both stochastic 
volatility dynamics and local volatility ones. This approach has to our knowledge never been studied 
in a stochastic interest rates framework but gives positive results for short dated options when interest 
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et al. 2007, Tavella et al. 2006] ). In the constant interest rates case, once the local volatility surface 



is available, the new mixed volatility can be computed by multiplying this local volatility with a ratio 
of integrals that depend on the joint density of the FX spot and the stochastic volatility. This density 
can be determined by numerically solving the associated two-dimensional Kolmogorov forward PDE. 



The study of the local volatility and its calibration in a three-factor model with local volatility 
can therefore be motivated by hedging arguments but is also considerably useful for the calibration of 
hybrid volatility models. In this paper we derive the local volatility function in a three-factor model 
with local volatility where we have three sources of randomness: the FX spot with a local volatility 
and the domestic/foreign interest rates. 

In a one-factor Gaussian model, the local volatility surface is generally built by using the Dupire's 
formula where partial derivatives of call options with respect to strikes and maturities are calculated 
by finite differences and where the real implied volatility surface is an interpolation of a finite set of 
market call prices. In a three-factor framework with local volatility, the expression of the local volatil- 
ity becomes more complicated as it also depends on a particularly complicated expectation where 
no closed form expression exists and it is not directly related to European call prices or other liquid 
products. Its calculation can obviously be done by numerical integration methods or Monte Carlo 
simulations. To enable realizations of the numerical integrations you need the forward probability 
distribution of the spot FX rate and the domestic and foreign interest rates up to maturity which is 
a solution of a three-dimensional forward PDE. 



An alternative approach is to calibrate the local volatility from stochastic volatility models by 
establishing links between local and stochastic volatility. Extracting the local volatility surface from a 
stochastic volatility model rather than by using the market implied volatility surface presents several 
advantages. First, the market implied volatility surface can in practice only be an interpolation of a 
finite set of available market data. As a consequence a local volatility surface built from an approxima- 
tive implied volatility surface is often unstable. Stochastic volatility models can be calibrated by using 
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fast algorithms like Fast Fourier Transforms (FFT) (see for example Ahlip, 2008, Andreasen 2006 



van Haastrecht et al. 2009 1) and the local volatility surface extracted from the calibrated stochastic 



volatility model is really smooth. In this paper we will present some mimicking properties that links 
the three-factor model with a local volatility to the same model with a stochastic volatility rather 
than a local volatility. These properties will allow us to obtain explicit expressions to construct the 
local volatility surface. 

Finally, we derive a link between the three-factor model with local volatility and a hybrid volatil- 
ity model where the volatility of the spot FX rate mixes a stochastic volatility with a local volatility. 
Knowing the local volatility function associated to the three-factor model with local volatility, we 
propose a calibration method for the local volatility in the four-factor hybrid volatility model. 

Our objective is to expose theoretical results about the local volatility function and its calibration. 
Numerical implementation of the model to the long-dated FX options market will be studied in a 
subsequent paper. 

This paper is organized as follows: We begin by defining the three-factor model with local volatility 
in section [2] Then, in section [3] we derive the local volatility expression for this model. First, we 
derive from the three-dimensional Fokker-Plank equation for the forward probability density function, 
a "simpler" one-dimensional PDE. This PDE is used in the derivation of the local volatility function 
by differentiating European call price expressions with respect to the strike and the maturity. Section 
[4] is devoted to the calibration of this local volatility function. In section 4.1, we present two numerical 



approaches based on respectively Monte Carlo simulations and PDE numerical resolution. In section 



4.2 we obtain a link between the local volatility function derived in a three- factor framework and the 



one coming from the simple one-factor Gaussian model. Next, in section 4.3, we derive a link between 
the three-factor model with a stochastic volatility for the spot FX rate and the one where the spot 
FX rate volatility is a local volatility. This link provides a relationship between our local volatility 
function and future instantaneous spot FX rate volatilities. Finally, in section[5j we study an extension 
of the three-factor model with local volatility. We derive a link between the local volatility function 
associated to the three-factor model and the local volatility function present in the four-factor hybrid 
volatility model and propose a calibration procedure for the local volatility function associated to this 
four- factor hybrid volatility model. Conclusions are presented in section [6j 

2 The three-factor pricing model with local volatility 

In this paper we consider the three-factor model where the spot FX rate volatility is a deterministic 
function of both the time and the FX spot itself. This function is known as 'local volatility'. In this 
model, the spot FX rate S is governed by the following dynamics 



dS(t) = (r d (t) 



r f {t))S(t)dt + a(t, 5(t))5(t)dW| ,iW (t), 



(1) 



where the volatility of the spot FX rate is a local volatility a(t, S(t)) and where domestic and foreign 



interest rates denoted by and rj respectively follow a Hull- White one factor Gaussian model Hull 
and White, 1993 defined by the Ornstein-Uhlenbeck processes 
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J dr d {t) = [9 d (t) - a d {t)r d {t)]dt + a^dW^ (t), (2) 
\ dr f (t) = [0 f (t) - a f (t)r f (t)-p fS <Tf(t)a(t,S(t))}dt + a f (t)dWf RN (t), (3) 

where 9 d (t),a d (t),a d (t),9f(t),af(t),<jf(t) are deterministic functions of time. Equations ([!]), Q and 
([3]) are expressed in the domestic risk- neutral measure (DRN). 

Foreign and domestic zero-coupon bonds defined by Pi(t,T) = E*[e - /* n(s)ds | F t ], i = d,f 
satisfy the evolution 

f dP d (t,T)=r d (t)P d (t,T)dt-a d (t,T)P d (t,T)dW d DRN (t), (4) 
\ dP f (t,T) = [r f (t) + p fS a f (t,T)a(t,S(t))]P f (t,T)dt - a f (t,T)P f (t,T)dWf RN (t), (5) 

where their volatility structures are given by 

f-T 

*i(t,T) = ai(t)J e-f* ai M du ds, i = d,f. 



More details about stochastic interest rates in the Hull- White one factor Gaussian model and 
explicit expression for the price of zero-coupon bonds under this framework can be found in |Brigo 



and Mercurio 2006]. 

Finally, in this paper we assume that (Wg RN (i), W d R (t), Wf* (t)) is a Brownian motion under 
the domestic risk-neutral measure Q d with correlation matrix 

1 PSd PSf* 
PSd 1 Pdf 
\PSf Pdf 1 

3 The local volatility function 

In this section we derive the expression of the local volatility function by using the same technique as 



Dupire |Dupire, 1994 and Derman and Kani Derman and Kani, 1994 which consists in differentiating 



European call price expressions with respect to the strike and the maturity. 
3.1 Forward PDE 

Consider a world where both the spot FX rate as well as the domestic and the foreign interest rates 
are stochastic. Under the assumption of absence of arbitrage opportunities, and working under the 
i-forward measure Qt (where the domestic zero-coupon bond is the numeraire), the present value 
V(S(0),r d (0),rf(0),t = 0) of a derivative that pays off V(S(t),r d (t),rf(t),t) at time t is given by 



V(S(0), r d (0), 77(0), t = 0) = P d (Q, t)E* [V(S(t),r d (t),r f (t),t)] 

= Pd(0,t) J J J V(x,y,z,t)(j) F (x,y,z,t)dxdydz 
5 



where cf)F(x,y, z,t) corresponds to the i-forward probability density. 



Following the same approach as in 'Equity Hybrid Derivatives' Overhaus et al. , 2006] (readapt- 
ing the method in the context of FX market), we can get the following Fokker-Plank equation for 

<f>F(x,y,z,t) 



+ 



dt + {rd{t) ~ m t]) + dx + dy 

d[(9 f (t) - a f (t) r f (t) - PfS <T f (t)a(t, S(t)))M 1 d 2 [a 2 (t, S(t))S 2 (t)0 F ] 1 gV|ft)^] 

dz 2 dx 2 2 dy 2 

ld 2 [a 2 (t)<f> F } 3 2 [a(i,S(t))S(t)a d (tW F ] d 2 [a(t, S(t))S(t)a f (t)p Sf cf> F ] &[j d (t)* f (t) P4 f<t>i 



2 dz 2 dxdy dxdz dydz 

Equation (j6| is a forward PDE since it is solved forward in time with the initial condition at time 
t = given by 4>f(x, y, z, t) = 5{x — xq, y — yo,z — zq), where 5 is the Dirac delta function and xq, yo 
and zq correspond to the values at time t = of the spot FX rate, the domestic and foreign interest 
rates respectively. This initial condition means that at time t = we are sure that the spot FX rate 
S(0) equals xo, the domestic interest rate r d (0) equals yo and the foreign interest rate rj(0) equals zq. 

From equation Q we will derive a "simpler" one-dimensional PDE only involving the dependence 
in x of the forward density function. This equation will be useful in the derivation of the local volatility 
expression in section [3^2] 

Let us denote by qp(x,z,t) the function defined as the integral of <(>f(x, y, z, t) over the whole 
range of y 

/+oo 
4> F {x,y,z,t)dy. (7) 
-00 

Making the realistic assumptions that lim cj>F(x,y, z,t) = and that partial derivatives of <f>p 

y— >±oo 

with respect to x, y and z tend also to zero when y tends to infinity, we obtain, after integrating ^ 
with respect to y, the following two dimensional PDE 

= ^+ J(rdt)-fd(O,t))cPFdy+^(J(r d (t)-r f (t))S(t)0 F dy^ 

d[(9 f (t) - a f (t) r f (t) - p fS CTf(t)a(t, S(t)))q F ] 1 d 2 [a 2 (t, S(t))S 2 (t)q F ] 



+ 



dz 2 dx 2 



1 d 2 aj(t)q F d 2 a(t, S(t))S(t)a f (t)p Sf q F 



2 dz 2 dxdz 
Finally, we define pf(x, t) as the integral of qF(x, z, t) over the whole range of z 

/+oo r+oo r+oo 

q F (x,z,t)dz= / <p F (x,y,z,t)dydz. (9) 

-oo J — oo J~oo 
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(6) 



Assuming that lim q F (x, z,t) = and that partial derivatives of q F with respect to x and z tend 

z— »±oo 

also to zero when z tends to the infinity, then integrating ^ with respect to z reduces the PDE ^ 
into the following one-dimensional PDE 



dpF 
dt 



+ 



(r d (t)- f d (0,t))<p F dydz + 



d_ 
dx 



(r d (t) -r f (t))S(t)<t> F dydz 



l d 2 [a 2 (t,S(t))S 2 (t)p F } 
2 dx 2 



(10) 



3.2 The local volatility derivation 

In this section we outline how to derive the expression of the local volatility function associated to 
the three-factor model with local volatility by differentiating the expression of a European call price 
(C(K, t)) with respect to its strike K and its maturity t. In order to simplify our calculations, we will 
work with the forward call price C(K, t) defined by 



+oo 



(x — K)(j)p(x, y, z, t)dxdydz. 



(11) 



K 



Differentiating equation (11) twice with respect to K, we obtain 



dC(K,t) 

dK 
d 2 C(K,t) 
dK 2 



e r r+oo 

J J Ik ~^ F ^ x ' y ' z '^ dxdydz = ~ EQ '[ 1 {5(i)>^}] 5 



4>f(K, y, z, t)dydz = p F (K, t). 



(12) 
(13) 



Differentiating equation (11) with respect to the maturity and using equation (10) leads to 



dC(K,t) 
dt 




+C °(S(t)-K) d4>F{x ^ Z ^ dxdydz 



ST [s[i) - K)( t dx 



+oo 



(S(t)-K){- 



K 



(r d (t) - f d (0,t))cj> F dydz 







I I (r d {t) - r f (t))S(t)</> F dydz) + 



\ d 2 [a\t,S{t))S\t)p F ] 
2 dx 2 



}dx. 



Integrating by parts several times and using equation ( 13 ) we get 
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dC ft ,1:) = fd(0,t)C(K,t) + J J J^ 00 [r d (t)K-r f (t)S(t)} ( l )F (x,y,z,t)dxdydz 

+ -(a(t,K)K) QR2 
= f d (0, t)C(K, t) + [(r d (t)if - r/(t)5(t))l {s(t)>x} ] 

+ > WtlJO iO.». ,14) 

This leads to the following expression for the local volatility surface in terms of the forward call 
prices C(K, t) 

aHt R) = ~ fd(0, t)C(K, t) - EQ* [{r d {t)K - r f (t) S(t))l {s{t)>K} } 

° [ ' ' 1 K 2 d 2 C(K,t) ' [ ' 

2^ dK 2 

The (partial) derivatives of the forward call price with respect to the maturity and the strike can 
be easily calculated 

adiK.t) _ aiig] _ 8C(K,t) i +U0]t)d(KJh (16) 



dt dt dt P d {0,t) 

d 2 C(K,t) ^[gggl _ 1 d 2 C(K,t) 
dK 2 dK 2 P d (0,t) dK 2 



(17) 



Substituting these expressions ( |16| ) and (17) into (15), we obtain the expression of the local volatil- 
ity cr 2 (t, K) in terms of call prices C(K, t) 



dC(K,t) 
JZf-i ts\ — dt 



P d (0, t)E* [(r d (t)K - r f (t) S(t))l {s{t)>K} 



2 A dK 2 

Unfortunately, this extension of the Dupire formula is not easily applicable for calibration over 
the market since there seems no immediate way to link the expectation term with European option 
prices or other liquid products. However, we present in section [4] four different methods to calibrate 
this local volatility function. 



Finally, we underline the fact that when assuming deterministic interest rates, equation (18) re- 
duces to the simple Dupire formula corresponding to the one factor Gaussian case. This formula can 
easily be derived by a similar reasoning as above but in a one factor framework. 

2 (f K) ggg + (r d (t) - r f (t))K 9 -^ + r f (t)C(K, t) 

2 A dK 2 
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Remark 1 The market often quotes options in terms of implied volatilities crimp instead of option 
prices. Consequently, it is more convenient to express the local volatility in terms of implied volatilities 
than option prices. As the implied volatility of an option with price C(K,T) is defined through the 
Black- Scholes formula (C mkt (K,T) = C BS (K,T,ai mp )), the derivatives of call prices in equation (19) 
can be computed through the chain rule, and this leads to the following equation (see [Wilmott, 20061) 



a 2 (t,K) 



9 ^JHR +2{r d (0) - r f (0))Kta ir , 



' imp 



[1 + Kd+Vi^^) 2 + KHo. 



imp y qj{2 



rap \2 



dK 



(20) 



Using the same approach, the local volatility expression (IS) can be written in terms of implied 
volatilities Uj, 



'imp; 



a 2 [tK) = a . e-^(°) f g(0){^(d + )(a, mp + 2t%g) - 2y/ir f (0)Af(d+)} + 2y/i{r d (0)Ke- r *WW(d-) + E} 
^ e-^/(°)*5(0VV'(d + ){(l + Kd + Vt^) 2 + KHa imp (%^ _ d+( ^E)2^ )} 

(21) 

where 



E = P d (0,t)B^[(r d (t)K - r f (t) S(t))l m)>K} ], 



d± 



log S ^ + (r d (0)-r f (0)± a ^)t 



a 



imp 



Vi 



M{x) 



x l -jfl 

e 2 dz, 

oo 



2vr 



1 - a 2 
N (x) = ^=e 2 . 
V2tt 

4 Calibrating the Local Volatility 

Before using a model to price any derivatives, it is usual to calibrate it on the vanilla market which 
means that you should be able to price vanilla options with your model such that the resulting im- 
plied volatilities match the market-quoted ones. More precisely you have to determine all parameters 
present in the different stochastic processes which define the model in such a way that all European 
option prices derived in the model are as consistent as possible with the corresponding market ones. 

The calibration procedure for the three-factor model with local volatility can be decomposed in 
three steps: (i) Parameters present in the Hull- White one- factor dynamics for the domestic and foreign 
interest rates, 6 d (t), a d (t), a d (t),6f(t), ctf(t), o-f(t), are chosen to match European swaption / cap-floors 
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values in their respective currencies. Methods for doing so are well developed in the literature (see 



for example IBrigo and Mercuric 2006 ). (ii) The three correlation coefficients of the model, psd,PSf 
and pdf are usually estimated from historical data, (iii) After these two steps, the calibration problem 
consists in finding the local volatility function of the spot FX rate which is consistent with an implied 



volatility surface. In |Piterbarg, 20061, Piterbarg derives an approximative formula for European call 



prices in the three-factor model where the local volatility function for the spot FX rate is a parametric 
function of the form a(t,S(t)) = ^(^)(rpy)^ _1 > where v{t) is the relative volatility function, (3(t) 
is a time-dependent elasticity of variance (CEV) and L{t) is a time-dependent scaling constant. The 
calibration procedure of Piterbarg consists in determining the functions v{t) and (3{t) such that when 
pricing any European call with his approximative call valuation formula, the price he gets is as close as 
possible to the market call price. Next sections will be devoted to other calibration methods for the lo- 
cal volatility function a(t, S(t)) based on the exact expression of this function in a three-factor context. 



4.1 Numerical approaches 
4.1.1 A Monte Carlo approach 



In this section we present a first calibration method for the local volatility expression (18) derived in 
section^ In this approach, the expectation, ~Efi T [(r d (T)K — rf(T)S(T))l{ S (T)>K} ] is approximated 
by using Monte Carlo simulations up to a fixed time t = T. To calculate numerically this expectation 
we have to simulate the FX spot rate S(t) and both the domestic and foreign interest rates up to time 
T starting from the initial market prices S(0), r d (0) and ry(0) respectively. Since the expectation is 
expressed under the measure Qt, we have to use the dynamics of S(t), r d (t) and rf{t) under this last 
measure, 



r dS(t) = [r d (t) - r f (t) - a(t, S(t))a d (t)b d (t, T)p Sd ]S(t)dt + a(t, S(t))S(t)dW£ F (t), 
I dr d {t) = [9 d {t) - a d r d (t) - a 2 d {t)b d (t,T)]dt + a d (t)dWj F (t), 

{ dr f {t) = \0 f (t) - a f r f (t) - p fS a f (t)b f (t, T)a(t, S(t)) - a d (t)b d (t, T)a f (t)b f (t, T)] Pdf dt + a f (t)dWf F (t). 

where bj(t, T) = ^-(1 — e~ a ^ T ~ t ' > ), j = d, f (assuming that ay are constant). 

The idea of the Monte Carlo method is to simulate n times (i.e. n scenarios) the stochastic variables 
S(t), r d (t) and rf(t) up to time T, by using for example Euler discretisations. The expectation is 
approximated by: 

E^[(r d (T)K - r f (T)S(T))l {S{T>K} ] =i - f>^(T)if - r)(T)S*(T))l mT)>K} (22) 



n . 
i=i 



where i corresponds to the i -scenario i = 1, n. 

As we have to know the local volatility function up to time T to simulate the path for S(t) and 
rf(t), the only way to work is forward in time. To begin, we have to determine the local volatility 
function at the first time step T = T\ for all strike K. At this first step we assume that the initial 



local volatility is equal to the deterministic local volatility given by equation (19). Note that this 
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local volatility is directly obtained by using market data (see equation (20)). More precisely, by this 
choice, we assume that for a "small time period", interest rates are constant and in this case, the 
local volatility expression (18) reduces to (19). Knowing that local volatility function we can simulate 
S(T\) and r(T\). Then we can compute the expectation E^ T i \(r d (Ti)K — rf(Ti)S(Ti))lrs(Ti)>K}] f° r 
all K by using: 



1 n 

[{r d (T x )K - r / (T 1 )5(T 1 ))l {5(Tl)> ^ } ] - - ^(r^K - r^S^T,))!^^ (23) 



n 
i=l 



This allows us to get the local volatility expression at time T±, a 2 {T\, K), for all strike K. 

Following the same procedure we can easily calibrate the local volatility at time T2 by using the 
local volatility obtained at time T\ and also the simulated path until time T\. Following this procedure 
we are able to generate the local volatility expression up to a final date T = 

4.1.2 A PDE approach 

The strategy is to solve the forward equation ^ forwards one step at a time, starting with a local 
volatility <j(0, S(0)) at time to = 0. At the first time step t\ = to + At, we can generate the forward 
joint transition densities 4>f(%, U, z, to + At) by solving the forward PDE using the initial condition 
(Pf(x, y, z, to) = 6(x — xo,y — yo,z — Zq) at time to = 0. Knowing the forward joint transition densities 
we can calculate the expectation in (18) namely E^* [(r d (t)K — rf(t) S(t))l{ S (t)>K} ] a * this point t\. 
This allows us to calculate from equation (18) the function (j(t\,S(t\)) at this time step. Following 



this procedure through time, we generate simultaneously both the expectation in (18) and the local 
volatility function <r(t, S(t)). At each time step we generate also all domestic forward joint transition 
densities <j>F(S(t),rd(t),rf(t),t) by solving the forward PDE ffify. 



4.2 Comparison between local volatility with and without stochastic interest rates 

In this section we present a way to calibrate the extension of the local volatility function derived in 
section [3] in a three- factor framework by using the local volatility showing up in the simple one- factor 
Gaussian model. 

Assuming deterministic interest rates, the FX spot follows the diffusion equation 

dS(t) = (r d (t) - r f (t))S(t)dt + a lf (t, S(t))S(t)dW° RN (t) (24) 
where the local volatility function denoted by a\f is given by equation (19 ^ 



j (t K) _ gig + K(f d (0, t) - f f (0, t)) + //(0, t)C(K, t) 

However, if we consider the three- factor model with stochastic interest rates, the local volatility 



function is given by equation ( 18 ) 



in the case of deterministic interest rates fd(0,t) = rd(t) and //(0,t) = rf(t) 
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ac{K,t) 



, .. , ^ ~ P*(0, t)¥fi'[(Kr d (t) - r f (t)S(t))l {s(t)>K} ] 
<W*) = lT ^C(K,t) • ( 26 ) 

2 A OK 2 

Using the fact that under the T-forward measure for T > t, we have: E^ T [r d (T) | Ft] = f d (t,T) 
and E Qr [r/(T) | T t ] = ff(t,T) we notice that 



E^[r d (t)l {s{t)>K} ] = E^[r d (t)}E^[l {s{t>K} ] + Cov^[r d (t), l { s(t)>K}} 

= fd(P, t)(- ^ 9C ^ K) ) + Cov^ [r d (t), l {s{ t)>K}} (27) 

where Cov^ 1 (X, Y) represents the covariance between two stochastic variables X and Y with dynamics 
expressed in the t-forward measure Qt- We also have 



E^[r f (t)S(t)l {s{t)>K }} 



B^[r f (t)(S(t) - K)l {s{t)>K} ]+B^[r f (t)Kl {s{t)>K} ] 
E Q *[r f (t)]E Qt [(S(t) - K) + ] + Cov Qt [r f (t), (S(t) - K)+] 
+K(E^[r f (t)]E^[l {s{t)>Ky ] + CovO*[r/(t), l{S(t)>K}]) 



+K(f f (0,t)(- 



dC(t,K), 



P d (0,t) dK 



+ Cov^[r f (t),l {S{t)>K} ]). 



(28) 



Substituting expressions (27) and (28) in equation (26), one finds the following interesting relation 
between the simple Dupire formula (25) and its extension (26) 



2 u 2 u KP d (0,t){C°v Qt lr f (t) - r d (t), l { s(t)>K}} + £Cav«* [»■/(*), (S(t) ~ 

2^ 

(29) 

Assuming that quantities Cov^' [rf(t), (S(t) — K) + ] and Cov^*[r^(t) — r d (t), l{s(t)>K}] are ex- 
tractable from the market, equation ( |29[ ) shows the corrections to make to the tractable Dupire local 
volatility surface in order to obtain the local volatility surface which takes into account the effects of 
both domestic and foreign stochastic interest rates. 



4.3 Calibrating the local volatility by mimicking stochastic volatility models 

In this section we consider the three factor model with a stochastic volatility for the spot FX rate and 
we show how to connect this model to the one where the spot FX rate volatility is a local volatility. 
More precisely, we sketch the derivation of an expectation relationship between local volatilities and 
future instantaneous spot FX rate volatilities. 
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Consider the following domestic risk neutral dynamics for the spot FX rate 



dS(t) = (r d (t) ~ r f (t))S(t)dt + 7 (t, u(t))S(t)dWg RN (t) (30) 

where v{t) is a stochastic variable which provides the stochastic perturbation for the spot FX rate 
volatility. Common choices for the function j(t, v(t)) are v(t), exp{^ v(t)) and sj v(t) . The stochastic 
variable v(t) is generally modelled by a Cox-Ingersoll-Ross (CIR) proce ss (mean-reve rsion in the drift 
and a volatility dependent of \J v(t) as for example the Heston model 



Heston 



1993 



) or in the first 

case by a Ornstein-Uhlenbeck process (OU) (mean-reversion in the drift and a volatility independent 



from u(t) as for example the Schobel and Zhu Schobel and Zhu, 1999 stochastic volatility model) 
which allows for negative values of v (t) . 

Applying Tanaka's formula to the convex but non- differentiable function e~ h rd ^ ds (S(t) - K)+ 
leads to 



ti rd W ds {S(t)-K)- 



r-t 

(S{0) - K) + - / r d {u)e-^ r ^ s)ds {S{u) - K) + du 
Jo 

+ f t e -I u r d (s)ds 1 {s{u>K}dSu + l f e-S>^ ds dLK{S) 
Jo 1 Jo 



where L^{S) is the local time of S. Since S is a continuous semimartingale, then almost-surely (see 
|Revuz and Yorj|2001| ) 



Lf{S) = lim - f l [KtK+t] (S(s))d <S,S> S . 

elO e Jq 



(31) 



Using the domestic risk neutral diffusion for the spot FX rate (30), one obtains 
e -/ofi(»)^(5(t)-if)+ 



= (S(0)-K)+ + K f r d (u)e-fords)d Sl du 
Jo 

Jo 

+ [\-^^ d n {s{u)>K}1 (u,u(u))S(u)dW^ RN (u) 
Jo 

+ / e -^^(s)ds dL K^ s y 
2 Jo 

Assuming that the function e~fo rd ^ ds l^ s ^ >K y^f(u,u(u))S(u) is a member of the class H 2 , 

namely the measurable and adapted functions / such that E^ d [/ f 2 (s)ds] < oo, we get when taking 
the domestic risk neutral expectation of each side 

t I' 1 
E Qd [e-W s ) rfs (S(t) -K) + \ = B Qd [(S(0)-K) + ] + K E Qd [r d (u)e~ founds, j du 

Jo 

- ['^[e-ford^ {u) s{u)]du 

Jo 
ft 1 

+ / -¥P d [e- ^ rd ^ ds dL^ {S)]. 
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Differentiating this equation leads to 



dC(K,t) = KB^[r d (t)e-tir d (s)ds 1{s{t)>K} ]dt _ E Q d[e - /* ^W^i {5(t)>K}r/(t) S (t)]dt 

+ -E Qd [e~ ti rd{s)ds dL? {S)). 
2 



Using characterization (31 ) for the local time with d < S, S >t= 7 {t, v{t)) S (t)dt, we obtain 



dC(K,t) = KE^[r d (t)e-fora(s)d Sl{s{t>K} ]dt -E^[e~ ^ r ^ ds l {s{t)>K] r } (t) S(t)]dt 

+ ^hmE^[il [ ^ +e] (5(t))e-/o^WV(t,Ki))^ 2 W]^- (32) 



Here, the last terms of the equation (48) can be rewritten as 



hm^E^[l [m] (5(t)) e -/o^(^V(t,Ki))^(i)] 

hm-E^[E^[ 7 2 (t,Kt))e-i'o^(^ | S(t)]l [K , K+e] (S(t))S 2 (t)] 



eio e 



E Qd [-f 2 (t,u{t))e-fo r ^)ds | S y) = K] Pd (K,t)K 2 
EQi[-f 2 (t,v{t))e-ti r d(s)ds | = d 2 C{K,t) 



EQd[ e -fo r ^)^ | 5(f) = K] 



dK 2 



K 2 . 



(33) 



Substituting equation (49) in (48), one obtains 



E Qd [-f 2 (t, v(t))e-ti r ^ ds | S(t) = K] 
BQd[e-fo r "^ ds | S(t) = K] 



at 



EQ*[e-fo r ^ d %Kr d (t)-r f (t)S(t))l {s(t)>K} ] 



1 K 2 9 2 C 
2""- dK 2 



(34) 



The right hand side of this equation equals the local volatility in the three-factor model where the 
expectation is expressed under the domestic risk-neutral measure. Therefore, if there exists a local 
volatility such that the one-dimensional probability distribution of the spot FX rate with the diffusion 
([!]) is the same as the one of the spot FX rate with dynamics (30) for every time t, then this local 
volatility function has to satisfy the following equation 



a 2 (t,K) 



EQ*[-f 2 (t,v(t))e-fo rd ^ ds | S(t) = K] 
EQd[e~fo r d( s ) ds | S(t) = K] 



(35) 



Stochastic volatility models belong to the class of incomplete market models, and hence, we have 
to precise that this condition holds if we assume that the domestic risk neutral probability measure 
Q d used in both stochastic and local volatility framework are the same. This equation extends the 
result obtained by Dupire [Dupire 2004| and Derman and Kani Derman and Kani 1998 to the case 
where domestic and foreign interest rates are stochastic. We notice that when we assume deterministic 



interest rates, equation (51) reduces exactly to their result 
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a\t,K) = VP*[ 1 \t,v{t))\S{t) = K]. 



(36) 



Finally, we reduce equation (51) to a conditional expectation under the t- forward measure 



a\t,K) = B^[ 7 \t,v(t)) \ S(t) = K]. 



(37) 



The local volatility o~ 2 (t,S(t) = K) is, therefore, the conditional expectation under the t-forward 
measure of the instantaneous spot FX rate volatility at the future time t, contingent on the spot FX 
rate level S(t) being equal to K. 



Equation (59) give us a new way to calibrate the local volatility function. A numerical approach 



is to begin with solving the forward PDE Q by using for example finite differences or finite elements 
methods, to determine the forward probability density <pF and afterwards using numerical integration 



methods in order to calculate the conditional expectation (59). 



Conditional expectations are difficult to compute by using traditional Monte Carlo simulations 
since paths generated by the simulation will miss the event involved in the conditional expectation. 
However, Malliavin integration by parts allows to obtain different representations of such conditional 
expectation that can be compute efficiently by Monte Carlo simulations (see |Fournie et al. 2001| ). 



Remark 2 If we assume independence between the spot FX rate and its volatility, the local volatility 
function is given by a 2 (T,K) = E^ T [7 2 (T, v{T))\ . In some particular cases it is possible to derive 
closed-form solutions for this expectation. Consider for example the three-factor model with local 
volatility where domestic and foreign interest rates have Hull and White dynamics (see section^ and 
we want to find the local volatility function by using the mimicking property ( 59 ) derived from a four- 



factor model with the same interest rates dynamics and with Schobel and Zhu dynamics for the FX 
spot volatility (v(f) ), 



dS(t) -- 

dr d {t) 

dr f {t) 

[ du(t) = 



r f (t))S(t)dt + v(t)S(t)dWg RN (t), 
-a d r d (t)]dt + a d dW^ RN {t), 



(r d (t) - 

= Mt) 

= [9 f (t) -a f r f (t) - PfS a f v{t)]dt + a f dWf RN (t), 
k[\-v(t)] dt + £dW? RN (t), 



Under the T -Forward measure, the Schobel and Zhu dynamics (41) become 



(38) 
(39) 
(40) 
(41) 



dv{t) = k{k{t) - v{t))dt + £ dWi * (t) 
where A(t) = A - b d {t,T) = i(l - e -«d(r-t)) 



(42) 



Integrating equation (42), we obtain, for each t <T 



v{T) = v(t)e~ k ( T ~^ + [ T kA(u)e- k( ~ T - u Uu+ f dW? F {u), 
Jt Jt 
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(43) 



so that v(T) conditional on Tt is normally distributed with mean and variance given respectively by 



E«>(T)|Ji] = u(t)e~ k ^ + (A - P -^A)(l - e-KT-t)) + P*>°dS (1 _ e _( ad+fe )(T- t)) 

Y^ T [v{T)\T t \ = £(l-e- 2fe ( T -*)) 
Finally, the closed form solution for the local volatility function is given by 



a 2 (T,K) = E Qt [v 2 (T)] = (E QT [^(T)]) 2 + Var QT [z/(T)] 

= f Ki)e" fcT + (A - ^*)(1 - e~ fcr ) + ^4^ (1 - e-^+^jV + £(1 - e" 2fcT ) 

(44) 



5 Hybrid volatility model 

In a simple local volatility model, the instantaneous volatility of the spot FX rate is a deterministic 
function of time and spot FX level. As a consequence, these models are suitable for pricing derivatives 
in situations where the spot FX volatility is strongly correlated to the spot FX market level itself. 
However, in the FX option's market, the volatility seems to exhibit some stochastic behavior especially 
in the long dated market. In this section we study an extension of the three- factor model with local 
volatility which incorporates stochastic behavior in the spot FX volatility without deleting the local 
volatility one. 



In this section we consider a hybrid volatility model where the volatility of the spot FX rate is a 
combination of a stochastic and a local volatility function. More precisely, the volatility for the spot 
FX rate corresponds to a local volatility crLOC2(t, S(t)) multiplied by a stochastic volatility j(t, v{t)) 
where v(t) is a stochastic variable. This gives the following four- factor model with local volatility 



' dS(t) = (r d (t) - r f (t))S(t)dt + o- LO C2(t, S(t)) 7 (t, v(t))S(t)dWl? RN (t), (45) 
dr d (t) = [6 d (t) - a d (t)r d (t)}dt + a d (t)dW d DRN (t), 

dr f (t) = [0 f (t) - a f (t)r f (t) - p fS (Tf(t)a L oc2(t, S(t)) 7 (t, v(t))]dt + a f (t)dWf RN (t), 
k du(t) = a{t, u(t))dt + 0(t, v(t))dW? RN {t). 

We will show how to connect this model to the one where the spot FX rate volatility is a pure 
local volatility. In this case the volatility of the spot FX rate is modelled by a local volatility denoted 
by o-LOCi(t,S(t)) 



( dS(t) = (r d (t)-r f (t))S(t)dt + a L oci(t,S(t))S(t)dWl? RN (t), 
| dr d (t) = [9 d (t) - a d (t)r d (t)]dt + a d (t)dW? RN (t), 

{ dr f (t) = [6 f (t) - a f (t)r f (t) - P fsaf(t)a LO oi(t, S(t))]dt + a f {t)dWf RN {t). (46) 

Consider the domestic risk neutral dynamics for the spot FX rate in the Hybrid volatility model 
(equation(45)) 
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dS(t) = (r d {t) - r f (t))S(t)dt + a LOC 2(t, S(t))j(t, v{t))S{t)dW^ HN (t) (47) 

where v{t) is a stochastic variable which provides the stochastic perturbation for the spot FX rate 
volatility. Applying Tanaka's formula to the non- differentiable function e~ fa r ^ ds (S(t) - K)+ and 



following the same steps as in section 4.3 leads to 



dC(K,t) = KB^[r d (t)e'fo r ^ ds l {s{t)>K} ]dt - [e~ fo r ^ ds l {s{t)>K} r f (t)S(t)]dt 

+±limE«'£l [Jr>JrH (S(t^ (48) 



Here, the last terms of the equation (48) can be rewritten as 



limiE^[l[^ +e] (S(^ 

= lim^E^[E^[ 7 2 (t,Kt))e-^^(^ | S(t)]l [K , K+£] (S(t))al OC2 (t, S(t))S 2 (t)] 
= EQ*tf(t,v(t))e-fi r *M d > | S(t) = K] Pd (K,t)a 2 LOC2 (t,K)K 2 

] -a 2 LOC2 (t,K)K\ ( 19) 



E^[7 2 (t,i/(t))e-/o* r <»W d - | S(t) = K] d 2 C{K,t) , 



EQd [e" £ | 5(f) = K] dR2 



Substituting equation (49) in (48), one obtains 



B^[ 7 2 (t,u(t))e-fo r ^ ds | S(t) = K] _ ^-BQ^[e-fo r ^ ds (Kr d (t)-r f (t)S(t))l {s(t)>K} ] 



The right hand side of this equation equals the local volatility in the three-factor model where the 
expectation is expressed under the domestic risk-neutral measure. Therefore, if there exists a local 
volatility such that the one- dimensional probability distribution of the spot FX rate with the diffusion 



(46) is the same as the one of the spot FX rate with dynamics (47) for every time t, then this local 



volatility function has to satisfy the following equation 

°ioc«,K) - ^.io^vT^TTT* 1 (51) 

EQd[ e -fo r ^ ds | S(t) = K] 
= a 2 LOC2 (t,K)B^[l 2 (t,u(t)) | S(t) = K\. (52) 

Hybrid volatility models belong also to the class of incomplete market models, and hence, we 
assume that the domestic risk neutral probability measure Q d used in both hybrid and pure local 



volatility framework are the same. This equation extends the result obtained by Madan Madan et al 



2007 to the case where domestic and foreign interest rates are stochastic. We notice that when we 



assume deterministic interest rates, equation (51) reduces exactly to their result 
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a 2 LOC1 (t,K) = a 2 LOC2 (t,K)E^ 2 (t,v(t)) | S(t) = K]. (53) 

The local volatility a LOC2 (t, K) is, therefore, given by the local volatility a 2 LOCl {t, K) divided by 
the conditional expectation under the i-forward measure of the instantaneous spot FX rate volatility 
at the future time t, contingent on the spot FX rate level S(t) being equal to K. 



cr LOC2( t , K ) 



EQttf(t,v(t))\S(t) = K] 



(54) 



Equation (54) give us a way to calibrate the local volatility function <JLOC2(t, K) knowing the local 
volatility function aLoci{t^)- A numerical approach is to begin with solving the forward PDE ^ 
by using for example finite differences or finite elements methods, to determine the forward probability 
density cpp an d afterwards using numerical integration methods in order to calculate the conditional 



expectation (54) (see subsection 5.1). As mentioned in section 4.3, an other approach is to use Malli- 
avin integration by parts to obtain different representations of the conditional expectation that can 
be computed efficiently by Monte Carlo simulations (see |Fournie et al. 2001 1 ) . 



Remark 3 // we assume independence between the spot FX rate and its volatility, the local volatility 
function is given by a 2 LOC2 {t,K) = ^^^i^~^yyj- Consider the hybrid model with Schobel and Zhu 
dynamics for the stochastic component of the FX spot volatility (v(t)), 



r f (t))S(t)dt + o-LOC2(t, S(t))v(t)S(t)dW° RN (t), 
a d r d (t)]dt + a d dW° RN (t), 



dS(t) = (r d (t) 

dr d {t) = [6 d (t) — u. d i d \b)\u,L -r u d uw d 
dr f (t) = [9 f (t) - a f r f (t) - p fS a f u{t)]dt + a f dWf RN (t), 



k dv{t) = k[X - u(t)} dt + £dW° RN {t) 



(55) 
(56) 
(57) 
(58) 



In this particular case we have a closed form solution for the local volatility function o~loC2, 
(knowing the local volatility function a loci) given by 



a LOC2( t , K ) 



(v{t)e~ kT + (A 



P 2 LOC\^ K ) 



a d k >\ 



-kT\ 



+ 



a d (a d +k)) 



(1 



-(a d 



-2kT\ 



.(59) 



In this section we obtained a link between the pure local volatility model (46) and the hybrid 



volatility model (45) such that the one-dimensional marginal distributions of St is the same for ev- 



ery t in both models. Using this result for the calibration of the local volatility o~LOC2(t, S(t)) (see 
section 5.1) will guarantee that prices of Vanilla options will be the same in both models since these 



products are fully determined by the conditional probability of the spot FX rate at the maturity time 
T conditional to the actual spot FX rate value. However, when we price exotic and path dependent 
options where payoff depends on intermediate spot values, pricing requires the product of all condi- 
tional probabilities. By using this link, we have never imposed that these probabilities are the same. 
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As a consequence, prices coming from the hybrid volatility model (45) are different from those coming 



from the pure local volatility one (46) as well as the pure stochastic volatility one (30) in the case 
of exotic derivatives. The advantage of the hybrid volatility model is that we have a combination of 
a local and a stochastic volatility effect in the volatility of the spot FX rate which is more realistic. 
Consequently, hybrid volatility models are generally more consistent with the exotic option's market 
after being calibrated with respect to the vanilla market (see Lipton and McGhee, 2002 1). 



5.1 Calibration 

In this subsection we are interested in the calibration of the local volatility (JLOC2(t, S(t)) of the four- 
factor hybrid volatility model knowing the local volatility (JLOCi{t-, S{t)) associated to the three-factor 



model that we have studied in the first part of this paper. In [Madan et al. , 2007 , the authors obtain 



a fast and exact calibration of vanilla options in the case of equity options with deterministic interest 
rates. Based on results obtained in this section, we extend this calibration procedure to the case of 
FX derivatives in the settings of stochastic interest rates. 

Remember that parameters present in the Hull- White one-factor dynamics for the domestic and 
foreign interest rates, 9d{t),a.d{t),ad{t).,6 f(t),af(t),a j(t), are chosen to match European swaption / 
cap-floors values in their respective currencies. The four correlation coefficients of the model, psd, PSf-, 
Pdf and psu are deduced from historical data. In this hybrid volatility model one also has to calibrate 
the process for v{t). The volatility of the variance i9(t,^(t)) can be deduced from historical time 
series with for example daily realized volatilities. If one uses a CIR process to model the stochastic 
variable i?(t, z;(t)), then it is of the form i!)(t,v(t)) = /3(t)y / z^(t), while if one models it by using an 
OU process, $(t, v{t)) is independent of v(t) and given by $(t,z/(t)) = fi{t). Usually, practitioners 
use a mean-reverting form for the drift of the variance, a(t, v{t)) = \{t) — Hi{t)v{i). One also has to 
calibrate the local volatility (TLOCi(t, S(t)) associated to the three-factor model with a local volatility 
only. This calibration can be realized using methods developed in section |4j Once all these parameters 
have been determined, one is able to calibrate the local volatility function crioc2(t, S(t)) associated 



to the hybrid volatility model by using equation (54), namely 



J2 i. ts\ _ u LOCI 

°LOC2(t,K) - E Q tb 2 {tiU{t))ls{t) = K y 



where the conditional expectation is by definition given by 



B^[ 1 2 (t,u(t))\S(t) = K] 

_ f™ 1 2 (t,v(t))MS(t) = K,r d (t),r f (t),v(t),t)dv 



(61) 



/ °° MS(t) = K, r d (t), r f (t), u(f), t)dv 
The t-forward joint density (f)F(S(t), rd(t), rf(t) } v[t), t) satisfies a four-dimensional Kolmogorov 



forward equation that can be derived using the same method as in section 3.1 



The strategy is to solve this forward equation forwards one step at a time, starting with a local 
volatility <7zx>C2(0, S(0)) equal to <tloci(0, <S(0))/7(0, ^(0)) at time to = 0, since the volatility z/(0) 
is known and we have E Q * [7(0, v(0))\S(0) = K) = 7(0,1/(0)). At the time step t\ = t + At we can 
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generate the forward joint transition densities 4>f(x, y, z, u, to + At) by solving the forward PDE using 
the initial condition 4>f( x i V, z, u, to) = 8(x — xo, y — yo, z — zq, v — vq) at time to = 0. Knowing the 
i-forward joint transition densities we can calculate each integral in equation (61) to determine the 
function E^^ii, v(ti))\S(ti) = K] at this point. This allows us to calculate from equation (60) 
the function <JLOC2(ti, S(ti)) at this time step. Following this procedure through time, we generate 
simultaneously both the conditional expectation E^' [y(t, v(t))\S(t) = K] and the local volatility 
function &LOC2(t, S(t)). At each time step we generate also all forward joint transition densities 
(j) F (S(t),r d (t),rf(t),u(t),t) by solving the forward PDE. 



6 Conclusion 

We have derived the local volatility expression of the spot FX rate in a stochastic interest rates frame- 
work. Therefore this model is very promising in pricing and hedging long-dated FX derivatives which 
are more and more traded in the FX option's market. We have proposed four different approaches 
for the calibration of this local volatility function. First we have proposed two numerical approaches 
based on Monte Carlo methods and numerical resolution of a PDE, respectively. The third one is 
based on the difference between the tractable local volatility surface that exists in the context of 
deterministic interest rates and our generalized one. The last method consists in generating the local 
volatility surface from a stochastic volatility. More precisely, we have obtained an explicit link between 
the local volatility and the stochastic volatility and this in a stochastic interest rates framework. This 
method has the advantage to give a smoother and more stable local volatility surface than when it is 
built from market available option prices. 

Afterwards, we have considered an extension of the previous model which allows the volatility of 
the spot FX rate to have local and stochastic behavior. Indeed, we have studied a hybrid volatility 
model, where the volatility of the spot FX rate is the product of a local volatility and a stochastic 
volatility. We have obtained results useful for the calibration of this new local volatility with respect 
to the stochastic volatility and the local volatility studied in the first part of the paper. 

The three-factor model with stochastic volatility and its extension will be useful in the pricing and 
risk management of long-dated FX derivatives for which it is especially important to consider the risk 
of both domestic/foreign interest rates and the risk linked to the FX spot volatility. Future studies 
will include numerical tests of the calibration methods on the FX market and comparison between 
their speed and appropriateness to fit the market implied volatility surface. The impact of stochastic 
interest rates and hedging performance of these models are also left for future research. 
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